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By Edmund E. Calleghen and John 8. Serafini 


SUMMARY 


The effects of existing frictional heating were analyzed to deter- 
mine the conditions under which ice formations on aircraft surfaces can 
be prevented. А method is presented for rapidly determining by means 
of charts the combination of.Mach number, altitude, and stream tempera- 
ture which will maintain an ice-free surface in ап icing cloud. The 
method can be applied to both subsonic апа supersonic flow. The charts 
presented are for Mach numbers up to 1.8 and pressure altitudes from sea 
level to 45,000 feet. 


INTRODUCTION 


As the operational speed of aircraft is increased through the tran- 
sonic region, the frictional heating available to prevent the formation 
Of ice on the aircr&ft becomes an important quantity. The operation of 
missiles апа interceptor aircraft is therefore possible under pre- 
determined values of altitude, stream temperature, and Mach number which 
will preclude icing. The set of flight conditions which provides а 
surface temperature of 520 Е for a particular point on a body traveling 
in an icing cloud is termed the “icing limit" for that point. The 
analysis presented in reference 1 relates the frictional heating, evapor- 
ation, and heat transfer for this particular condition. This analysis 
was applied to a symmetrical diamond profile airfoil at zero angle of 
attack and it was found that the most critical region of the airfoil in 
terms of icing limit was on the rear surface behind the midchord 
(shoulder). Although this region is not subject to direct impingement 
of water droplets, the water film caused by droplets impinging on the 
front surface runs back over the shoulder and extends the hazard of 
icing over a large percentage of the airfoil. Icing might therefore 
occur at any point on a body which is subject to droplet impingement or 


runback. 


Experimental verification of the calculations presented in refer- 
ence 1 was obtained at the NACA Lewis laboratory on a diamond airfoil of 
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5.50 inch chord and 6.82 percent thickness for a Mach number of 1.56 and 
several pressure altitudes. The good agreement between experimental and 
analytical results provides credence for the analytical method presented 
therein. 


The procedure used for calculating the icing limit in reference 1 
can be applied generally to obtain the icing limit-for any body provided 
the flow conditions about the body are known. The charts presented 
herein are based on the analysis of reference 1 for calculating the icing 
limit at any point on a body. These charts are similar to that shown in 
figure 1 of reference 1 but cover an altitude range from sea level to 
45,000 feet. Several auxiliary charts are also presented from which the 
pressure coefficient, velocity ratio, or local Mach number at a particu- 
lar point on the body-can be converted to the parameters utilized in the 


icing limit charts. 


2 ` SYMBOLS 


The following symbols are used in this report: 


Ср pressure coefficient, Еа T 
1/2 АД 
ср specific heat of air at constant pressure, Btu/(1b) (°F) 
e vapor pressure, lb/sq ft | 
k. coefficient of evaporation 
E coefficient of heat transfer м - 205 0 Шақ к 
L latent heat of vaporization, Btu/lb 
М Mach number _ _ 
Ша, molecular weight of air 
Mo molecular weight of water vapor 
p static pressure, lb/sq ft 
r recovery factor 
T static temperature, ов 


То, е minimum free-stream static temperature corresponding to ice-free 
4 condition on surface as defined in eq. (2), OR 
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V velocity, ft/sec 

Y ratio of specific heats of &ir, 1.400 

p density of air, slugs /cu ft 

Subscripts: 

О free-stream static conditions 

1 local conditions at edge of boundary layer 
S surface 


ANALYSIS 


The generalized heat balance in ап icing cloud for an insulated 
body heated only by frictional effects is given as follows: 


(1) Heat due to the frictional, or viscous, effects 
plus 
(2) Heat due to the kinetic energy of the water droplets 
plus 
(3) Heat of fusion 
equal 
(4) Heat lost by convection 
plus 
(5) Heat for evaporation of water 
plus 
(6) Heat required to raise temperature of water droplets from stream 
temperature to surface temperature. 
Hardy's relation (ref. 2) is obtained by equating the frictional 
term (1) to the sum of the convective term (4) and the evaporative term 


(5). The heat of fusion term (3) is equal to zero if an ice-free 
surface is defined as being fully wet at 32° F with no ice particles in 


& NACA TN 2914 


the surface water film. In addition, terms (2) and (6) are of the same 
order of magnitude and tend to offset each other for calculations of the 
icing limit. Both increase in value with increasing Mach number; large 
values of either term are found only at large values of liquid-water 
content and flight speed. Such a combination is highly unlikely in 
actual Plight since high speeds usually occur at high altitudes, and 
high altitudes imply small liquid-water contents with the exception of 
flight in cumulus clouds of great vertical extent. 


The results of reference 1 indicate that the relation of Hardy 
(ref. 2) applies in particular to a surface barely wetted by a very thin 
film, but is probably nearly correct for the whole range of speeds and 
liquid-water contents of practical interest for high-speed flight. 


The relation of Hardy is given in reference 1 as 


Kama L 
ی‎ = т (1 + 20220) ма 4.25 на Е (1) 


Because interest in this problem is fixed on a definition of flight cir- 
cumstances that provides local surface temperatures of 32 F, the terms 
Та» L, and eg are constants. Hardy observed that for the range of 


temperatures near 32° F the ratio of the evaporation coefficient k, to 
the heat-transfer coefficient k, is very nearly unity. Іп addition, 


the air in a cloud is fully saturated at the static or free-stream condi- 
tions and the free-stream vapor pressure ес is therefore equal to the 


saturated vapor pressure at the free-stream static temperature. If it 

is assumed that tbe flow about the body, outside the boundary layer, is 
accomplished with no change in phase, that is, no condensation or evap- 
oration, then Dalton's law o? partial pressures applies and 


Equation (1) сап therefore be rewritten in the form 


r 2 
Ті (1 + T × м?) 12.75 ео 
492 -------------- Ту < - 2776.5 - ----- 
To c ? (р. /po) Po ~ 12.75 Bo = ес 
(2) 


which is quite convenient for purposes ОҒ calculation because 
(т то, а) (2 + r= т м2) апа ру /Ро сап be readily calculated at any 
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point on a body for shock-free flow from the pressure or velocity dis- 
tribution, the stream Mach number, and the recovery factor r. When the 
Plow about the body contains a shock, then the parameters 


(11/19, с) (1 кг T м2) апа p1/Po can be obtained. from the local 


Mach number, the stream Mach number, the pressure coefficient, and the 
recovery factor г. 


PREPARATION OF ICING LIMIT AND AUXILIARY CHARTS 


Because eo is a function of То о, the solution of equation (2) 


must be accomplished by trial and error. In an effort to reduce the 
laboriousness of such calculations a number of charts have been made and 
are presented in figure 1. For each pressure altitude considered, the 
free-stream static temperature corresponding to an ice-free condition 


я -l 
То,с is plotted as a function of the parameter (11/10, 4 (a + L r ма 
by means ОР equation (2) for constant values of P,/Po- For each set of 


values of P/Po and (Ито, о) (1 + r= г ма) ‚ а free-stream static 


temperature exists which corresponds to a surface temperature of 32° Е 
and therefore represents the minimum free-stream static temperature for 
an ice-free surface Toc: The lower limit of the free-stream static 


temperature was considered to be -40° F ‚ Since it has been shown (refs. 5 
and 4) that supercooled water droplets are not likely to exist below this 


temperature. 


The ten icing limit charts presented (figs. 1(а) to 1(3)) were cal- 
culated for pressure altitudes from sea level (fig. 1(а)) to 45,000 feet 
(fig. 1(j)) in 5000 foot increments. The use of the icing limit charts 
requires the knowledge of the parameters (ғ,/т, қ (2 + = r м, 2) ала 

2 
ра /Ро: For shock-free flow both these parameters can be readily calcu- 
lated from the pressure coefficient Ср ог the ratio of local to free- 
stream velocity үу Мо and the stream Mach number. For most subsonic 


airfoils either the experimental pressure coefficient or the calculated 
velocity ratio is known. A variety of methods is available in the lit- 
erature for determining local values of Ср or Vi о for almost any 


body їп shock-free flow (refs. 5 to 9). 


Tn the transonic-supersonic regime, wherein a shock normally exists 
in the flow field, the flow field about a body is not defined by the 
pressure coefficient and stream Mach number alone and an additional 
parameter such as shock wave angle or Local Mach number is required. 
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Most of the results in the literature are presented in terms of local 
pressure coefficient and local Mach number for constant values of free- 
stream Mach number (ref. 10). 


The pressure coefficient may be related to the pressure ratio and 


the stream Mach number by the compressible flow relations for a perfect 
gas without the assumption of isentropic flow as follows: 


Рі - Ро  2(21 - ро) 2 [PL 
2 Povo T Po Mo To 
(3) 
m 2 ú г 
+ Mo Cy = 1 + 0.7 Mo“ Cp 


Hence this result applies to both shocked апа shock-free flows. 


The velocity ratio may be related to the pressure ratio and the 
stream Mach number in the following manner for shock-free flow: 


2 2 
E Ра а Ро Po Рі M к. 8)" ( a] 
V ~ 2 0, 
Ë PoYo za о Po My 
The static pressure - Mach number relation for isentropic flow is given 
by 


T Q bue 
EI 
1+ حي‎ №? 
Solving for M^ gives 
у-1 
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Therefore 


(4) 


y-1 2 ; 
The parameter NA d + — r M; ) is related to P,/Po» the 
stream Mach number, and the recovery factor in the following way for 
shock-free flow: 


т-1 
T 7 р T 
— ) + چ‎ = му) = (55) ( + be њу 
О,е 0 
but 
_ Y 
1.2 fF у 1.2 
Y-i = [ — r- 


Multiplying both sides by r, adding unity to each side, and rearranging 
terms yield 


Е T 
-l ,, 2 1 -1 2 
1 + r ем = (= ( + та но) = + (1 - m 


I p 
l Y-1 а\ _ 7-1 „2 =) 
T 26) (+ two) =+ G - 2) [E 
0.286 


= (1 + 0.2 Mo”) r + (1 - г) (ж) (5) 
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For shocked flows the pressure ratio can be obtained by equation (3), 


but the parameter (о, (2 + Іі r м2) сап по longer be related to 


the pressure ratio and stream Mach number as for shock-free flows 


(eq. (5)). 


The parameter (m /то,е) H + r= r м2) сап, however, ре related 
quite simply to the local and free-stream Mach numbers. Тһе temperature 


ratio T, /T. is given in terms of Mach number by the energy relations 
which apply to both shocked and Shock-free flows as follows: 
1 + X 
п 1 +Z мо 
To el. д 
therefore „л .. .. Я _ 
1-2 -1 
T L + Mo 
Т — с + r= r w?) = —— ——À 5 с + r= r ња) (6) š 
0,c 1+ г. М 


The pressure ratio ру/ро can be readily obtained by use of equa- 
tion (3) if. Му апа Cy are known. Тһе term мог Ср сап be calculated 


from experimental pressure distributions which are usually presented for 
each value of Ма in terms of Ср. 


Theoretical velocity distributions Рог many bodies are given in terms 
of the velocity ratio or the square of the velocity ratio. The conversion 
from velocity ratio to pressure ratio (eq. (4)) is shown in figure 2, 
where the pressure ratio is plotted as a function of velocity ratio for 
constant values of the stream Mach number. 


Figure 3 ed the conversion from pressure ratio to the parameter 


(ео, (2 T inci = г мо) as E by equation (5) for shock-free flows. 


The parameter | (2 + == т = x мо?) is plotted as a function of the 
stream Mach number Mo for абы ene values of the pressure ratio P- /Po' 
Figures 3(a), 3(b), апа 5(с) are given for values of the recovery factor 
r of 0.80, 0.85, and 0.90, respectively. 

The relation given by equation (6) is shown plotteá in figure 4 where 
the parameter (т /то,е) 1 me r м2) ів plotted ав а function of the 
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stream Mach number for constant values of the local to free-stream Mach 
number ratio М1/М0. Figures 4(a), 4(b), and 4(c) are given for values 


of the recovery factor of 0.80, 0.85, and 0.90, respectively. 


Discussion and Illustrative Examples 


Most aircraft components which are exposed to icing are of airfoil 
or streamlined shape. For this general class of bodies most of the 
cloud water droplets impinge near the leading edge or nose region. The 
icing hazard, however, is not limited to the region of direct impinge- 
ment, since if the impinged cloud droplets do not freeze (because of 
frictional heating) they create a surface water film which flows back 
from the impingement area. The whole region wetted by the film must 
therefore be considered as a region of possible icing hazard. 


The critical point on an airfoil or similar body (1.е., the point 
where the largest value of То,с is required for a particular value of 


stream Mach number) occurs at the minimum pressure point because the 
decreased pressure enhances the evaporation and results in a reduced 
surface temperature. If impingement occurs only near the leading edge 
or nose region, the evaporation in the region ahead of the point of 
minimum pressure may be sufficient to remove all the water film and, 

in such a case, the results obtained from charts l(a) to 1(j) would not 
apply. The actual critical point might therefore occur forward of the 
minimum pressure point; therefore, values of To ,с based on minimum pres- 


sure at а particular Mach number would be greater than actually required. 


The use of charts l(a) to 1(j) сап be shown most readily by means of 
several examples. The following example illustrates typical results in 
the subsonic speed range. 


Example I. - Calculation is made of the free-stream static temper- 
ature required for an ice-free surface as a function of Mach number for 
ап NACA 65-206 airfoil at 15,000 feet pressure altitude and jJ. amgle of 
attack. 


The velocity ratio Vi/Vo can be found for this airfoil by the 


methods and results presented in reference 5. The maximum velocity ratio 
(minimum pressure ratio) for this airfoil was determined to be 
VNo = 1.139 anā occurs at the 45-percent chord station on the upper 


surface. Assuming a value of the recovery factor of 0.85 and that the 
surface is wet at this point determines То,е in the following way: 
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Free-stream geri кора ыы Parameter Minimum free-stream 
Mach number ratio, р. /р T " static temperature 

My у dur ci" 1 1 + = ү M^ T у 

(fig. 2(b) апё То, о 2 В 0262 

V ) 21 
1/0 |818. 5(v), Po? Moj [fig. 1(a)] 

0.400 486.8 

. 500 485.8 

· 600 480.0 

‚ 700 475.6 


The v&lues of Toye obviously decrease with increasing Mach number as 


would be expected. The results of example 1 show, however, that even 
at а Mach number of 0.7, protection is provided at stream temperatures 


greater. than 169 Е. 


Example 2. - This example is presented to show typical calculations 
and results in the transonic speed range. Determination is made of the 
relation between stream Mach number and stream static temperature which 
will provide an ice-free surface at the midchord of an 8.8 percent thick 
circular arc airfoil at zero angle of attack for altitudes of 10,000, 
25,000, and 40,000 feet. The recovery factor г is 0.90. 


The values of pressure coefficient Co and local Mach number M 


can be obtained from reference 10 for stream Mach numbers of 0.848 to 
1.500 and are listed as follows: 


TO,c; OR 


Altitude, ft 


‚| (fig. 
1(c))| 1(f)) 
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The results of example 2 аге shown plotted in figure 5, where the 
stream static temperature corresponding to an ice-free condition is 
plotted as a function of the flight Mach number for the three altitudes 
considered. The effect of altitude is clearly illustrated. The free- 
stream static temperatures required for protection are higher at the 
higher pressure altitudes, as would be expected; since, if other condi- 
tions are equal, the decreased pressure enhances the evaporation rate. 


Lewis Flight Propulsion Laboratory 
National Advisory Committee for Aeronautics 
Cleveland, Ohio, December 4, 1952 
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Figure 5. - Variation of minimum free-stream static temperature for ice-free surface 
as function of free-stream Mach number for station at 50 percent chord of 8.8 percent 
thick circuler are airfoil. 
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